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We investigate canonical, phantom and quintom models, with the various fields being non- 
minimally coupled to gravity, in the framework of holographic dark energy. We classify them and 
we discuss their cosmological implications. In particular, we examine the present value of the dark 
energy equation-of-state parameter and the crossing through the phantom divide, and we extract 
the conditions for a future cosmological singularity. The combined scenarios are in agreement with 
observations and reveal interesting cosmological behaviors. 
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I. INTRODUCTION 

Nowadays it is strongly believed that the universe is 
experiencing an accelerated expansion, and this is sup- 
ported by many cosmological observations, such as SNe 
la [1], WMAP [2], SDSS [3] and X-ray [4]. These ob- 
servations suggest hat the universe is dominated by dark 
energy with negative pressure, which provides the dy- 
namical mechanism of the accelerating expansion of the 
universe. Although the nature and origin of dark energy 
could perhaps understood by a fundamental underlying 
theory unknown up to now, physicists can still propose 
some paradigms to describe it. In this direction we can 
consider theories of modified gravity [5], or field models 
of dark energy. The field models that have been discussed 
widely in the literature consider a cosmological constant 
[6] , a canonical scalar field (quintessence) [7] , a phantom 
field, that is a scalar field with a negative sign of the ki- 
netic term [8, 9], or the combination of quintessence and 
phantom in a unified model named quintom [10]. The 
quintom paradigm intends to describe the crossing of the 
dark-energy equation-of-state parameter wa through the 
phantom divide —1 [11], since in quintessence and phan- 
tom models the perturbations could be unstable as w\ 
approaches it [12]. 

In addition, many theoretical studies are devoted to 
understand and shed light on dark energy, within the 
string theory framework. The Kachru-Kallosh-Linde- 
Trivedi model [13] is a typical example, which tries to 
construct metastable de Sitter vacua in the light of type 
IIB string theory. Despite the lack of a quantum the- 
ory of gravity, we can still make some attempts to probe 
the nature of dark energy according to some principles of 
quantum gravity. An interesting attempt in this direction 
is the so-called "holographic dark energy" proposal [14- 
17]. Such a paradigm has been constructed in the light of 
holographic principle of quantum gravity [18], and thus it 
presents some interesting features of an underlying the- 



ory of dark energy. Furthermore, it may simultaneously 
provide a solution to the coincidence problem, i.e why 
matter and dark energy densities are comparable today 
although they obey completely different equations of mo- 
tion [16]. The holographic dark energy model has been 
extended to include the spatial curvature contribution 
[19] and it has been generalized in the braneworld frame- 
work [20]. Lastly, it has been tested and constrained by 
various astronomical observations [21]. 

In the present work we are interested in investigating 
various field models of dark energy, where the fields are 
non-minimally coupled to gravity [22, 23]. Such models 
have been shown to present significant features and we 
study them in the framework of holographic dark energy. 
In particular, we examine the current value of w\ and 
the realization of a recent crossing through the phantom 
divide —1 from above. Additionally, we investigate the 
possibility of a future WA-divergence [9, 24], and the spec- 
ification of the time that is it going to happen. The plan 
of the work is as follows: In section II we construct the 
cosmological scenarios of non-minimally coupled canon- 
ical, phantom and quintom fields, in the framework of 
holographic dark energy. In section III we examine their 
behavior and we discuss their cosmological implications. 
Finally, in section IV we summarize our results. 



II. NON-MINIMALLY COUPLED FIELDS IN 
THE FRAMEWORK OF HOLOGRAPHIC DARK 
ENERGY 



Let us describe briefly the holographic dark energy 
framework [14-17]. In this dark-energy model one de- 
termines an appropriate quantity to serve as an infrared 
cut-off for the theory, and imposes the constraint that the 
total vacuum energy in the corresponding maximum vol- 
ume must not be greater than the mass of a black hole of 
the same size. By saturating the inequality one identifies 
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the acquired vacuum energy as holographic dark energy: 

3c 2 



Pa 
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(1) 



with L the IR cut-off and c a constant which can be set 
to 1. Although the choice of L has raised a discussion in 
the literature [16, 25], in this work we will use the Hub- 
ble scale. Note that the aforementioned choice for the IR 
cut-off has been found to have problems in conventional, 
minimally-coupled frameworks [15], but this is not any- 
more the case if one considers non-minimal coupling, as 
we do in the present work. Finally, we mention that the 
extension of holographic dark energy in the presence of 
non-minimally coupled fields could possibly raise some 
theoretical questions, but we assume that such an exten- 
sion is valid. The detailed examination of this subject is 
left for a future work. 

In the following, we are going to investigate holo- 
graphic dark energy in the presence of canonical, phan- 
tom, or both fields, non-minimally coupled to gravity. 
The space-time geometry will for simplicity be a flat 
Robert son- Walker: 

ds 2 = -dt 2 + a(t) 2 (dr 2 +r 2 <m 2 ), (2) 

with a(t) the scale factor. 



A. Canonical field 

We first consider a canonical scalar field with a non- 
minimal coupling. This case has been partially investi- 
gated in [26], and here we extend it. The action of the 
universe is 



1 
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R 
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where k 2 is a gravitational constant. In the action we 
have added a canonical scalar field </>, which in non- 
minimally coupled to the curvature with coupling param- 
eter ^0. Although we could include a specific potential 
(quadratic or exponential), for the scope of the present 
work and for simplicity we keep the form (3) since our 
results can be easily generalized to these potential-cases. 
Lastly, the term Cm accounts for the matter content of 
the universe. 

The presence of the non-minimal coupling leads to the 
effective Newton's constant: 
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k 2 (i-&kV) 



(4) 



The Friedmann equations and the evolution equation for 
the scalar field are [22] : 

K 2 (p M +PA + \ft + 6^##) 



+ + 6£ (H + 2ff 2 ) = (6) 

PM + PA + 3.ff (p M + PA+PM+ PA) = 0, (7) 

where H — a/ a is the Hubble parameter. In these ex- 
pressions, pm and Pm are respectively the pressure and 
density of the matter content of the universe. Finally, 
PA and pa are the corresponding components of dark 
energy, which as usual is attributed to the scalar field. 
Since we use the Hubble scale to define holographic dark 
energy, that is we take L — if -1 , (1) can be written as 
Pa = 3(87rG e //) _1 iJ 2 , which, due to the effective nature 
of the Newton's constant (4), leads to: 



PA = -z (1 - ^V) H 2 . 



(8) 



We are interested in extracting power-law solutions of 
the cosmological model (5)- (7), in the case of a dark- 
energy dominated universe (pm,Pm -C 1). Thus, we are 
looking for solutions of the form: 

a(t) — a n t r 

m = <t>ot s *. (9) 

Insertion of these ansatzes in equations (5), (6) yields: 
1) + 3r.s + 6r(2r - 1)£ = 

(10) 



s + 12^r = 0. 



As we can easily see, the case of conformal coupling (^ = 
1/6) is not interesting since it leads to the trivial case 
r = s<f, = 0. Thus, for ^0^1/6 we obtain: 



1 



r = 



4 - 24& 
3£0 



(11) 



leading to: 



a(t) = a i 4 - 24 «* 

3«J, 



= 0o * ^ ■ (12) 
We can use expression (8) in order to acquire p\(t): 

PA(t) = -^r a (t- 9 -^K a ^- 9 ). (13) 
Substitution into (7) then straightforwardly provides pa- 

Pa(*) = \r [t- 2 (2 - 3r) + &K 2 $t 2 "*- 2 (3r + 2s - 2)] . 

(14) 

In expressions (13) and (14), r and s<f, are given by (11). 
Hence, we can calculate the dark energy equation-of-state 
parameter w\(t) as: 



w A (t) 



PA(t) = 5 
p A (t) 3 



\ . (15) 



t- 2s *-^m 
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Relation (15) allows us to determine both the of w\- 
cvolution, as well as its current value waq- I n order to 
express it in a more convenient form for comparison with 
observations, we can set the current values do = 1 and 
t = 1, and use rlni = In a = — ln(l + z) with z the 
rcdshift. Therefore, we acquire: 



wa(z) = --■ 



2 + 



e -24£,ln(l+z)_£ 0/t 202 



(16) 



Expression (16) provides w\(z) in terms of the coupling 
parameter £^ and the amplitude <pQ. We discuss the 
cosmological implications for various sub-classes of the 
present model, in section III. 



B. Phantom field 

In this subsection we consider a phantom field with a 
non-minimal coupling, that is a field with an opposite 
sign in the kinetic term in the Lagrangian [8, 9]. Such 
models are widely used in order to acquire u>a less than 
— 1. The action of the universe is 



^R - \ i a a 2 R + \g^d^d v a + C M 

(17) 

and the presence of the non-minimal coupling leads to 
the effective Newton's constant: 



8nG eff =k 2 (1-^kV) 1 . 



(18) 



The cosmological equations and the evolution equation 
for the phantom field are [8] : 



2 k 2 (pm + Pa - \o 2 + 6^ a Ha&) 
ti 



3(l-^K 2 a 2 ) 



= 



a + 3H& - 6£ CT (h + 2H 2 ^j a = 



(19) 



(20) 



PM + PA + 3H (p M + PA + PM + PA) = 0. (21) 

Similarly to the previous subsection, the use of the Hub- 
ble scale in the definition of holographic dark energy, and 
the effective nature of the Newton's constant (18), lead 
to: 



PA = ^(l-t^ 2 <? 2 )H 2 . 



(22) 



We examine power-law solutions of equations (19)- 
(21), in the case of a dark-energy dominated universe 
(pm,Pm «!)• Thus, we impose: 



a(t) = a t r 
a(t) =a t s ". 



(23) 



Insertion in equations (19), (20) yields: 

s a {s a - 1) + 3rs a - 6r(2r - 1)& = 
s a - 12£ ff r = 0. 



(24) 



As we can easily see, the case £ CT = —1/6 leads to the 
trivial case r = s a = 0. Thus, for £ CT ^ —1/6 we obtain: 



1 



•V = 



4 + 24£ CT 



1 + 6^' 



(25) 



leading to: 



a(t) = a t«+24« 



a(t) = cr 
Using (22) we acquire: 
3 



PA(t) 

and thus (21) gives: 



(t- 2 ^n 2 a 2 t 2 ^- 2 ) 



(26) 



(27) 



1 



PA(t) = —r [t- 2 (2 - 3r) + ^K 2 a 2 t 2s °- 2 (3r + 2s a - 2)] , 

(28) 

where r and s a are given by (25). We can calculate the 
dark energy equation-of-state parameter WA(t) as: 



PA(t) 
PA(t) 
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Finally, similarly to the previous subsection, we can ex- 
press (29) in terms of the redshift z obtaining: 



wa(z) = - + 8£a 



f 2 2 



e 24^ ln(l+z) _ ^ K 2 a 2 



(30) 



Relation (30) provides wa{z) in terms of the coupling 
parameter £ CT and the amplitude (To. We examine it in 
detail in section III. 



C. Quintom model 

In this subsection we consider the quintom cosmolog- 
ical scenario [10], that is we consider simultaneously a 
canonical and a phantom field, both with non-minimally 
coupling. As we have stated in the introduction, this 
combined cosmological paradigm has been shown to be 
capable to describe the crossing of the phantom divide 
wa = —1. The action of the model is: 



/ 



2k 2 2 



2 R-\^ 2 R- 



\gr%WA + \^"d^d v a + Cm 



(31) 
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and the presence of the non-minimal coupling leads to 
the effective Newton's constant: 

8nG eff = n 2 [1 - n 2 {^ 2 +£ t( j 2 )]~ 1 . (32) 

The cosmological equations and the evolution equation 
for the canonical and phantom fields are [10]: 



leading to: 



a(t) = a 



<t>(t)=(ht 1 " 6 ^ 

a(t) = c q t . 



(41) 



H — 



PM + PA + \4> 2 - \<J 2 + 6^/f# + 6&Haa 



3[l-« 2 (^ 2 +i> 2 )] 
+ 3U0 + 6£ (A + 2ff 2 ) = 



(33) 
(34) 



Note also that we could cquivalently express solutions 
(40) in terms of Finally, as expected, the choice ^ = 
1/6 gives —1/6 < £ CT , which is also non-physical. 
Using (37) we obtain: 



PA(t) 



(r 2 -&K 2 $t 2 »*- 2 -&KVgt 



1 J2,2s„-1\ 



g + 3^(7 - 6£ CT (# + 2ff 2 ) cr = (35) and thus (36) gives: 



(42) 



PM + PA + 3H (p M + PA + PM + PA) = 0. 



(36) 



As usual, the use of the Hubble scale in the definition of 
holographic dark energy, and the effective nature of the 
Newton's constant (32), lead to: 



PA = — [l-«: 2 (£^ 2 +^> 2 )] H 2 . 



(37) 



We examine power-law solutions of equations (33)- 
(36), in the case of a dark-energy dominated universe 
(Pm,Pm «!)• Thus, we impose: 



a(t) = a t r 

m = u s * 

a(t)=a t s °. 



(38) 



Substituting into (33), (34), (35) and requiring a solution 
for all times we get: 

s (s - 1) + 3rs + 6r(2r - 1)& = 
s*(s« - 1) + 3rs a - 6r(2r - 1)6, = 
s + 12&r = 

s a - 12^r = 0. (39) 

In the present quintom scenario with both fields non- 
minimally coupled, it becomes clear that the existence 
of non-trivial solutions requires a relation between the 
couplings £0 and Thus, for the physically interesting 
case £0 ^ 1/6, we obtain: 



1 



r = 



3& 



•V = 



(40) 



Pa(*) = — r [tr 2 (2 - 3r) + £ K 2 2 t 2s ^ 2 (3r + 2s - 2) + 
H*K 2 <7 2 t 2s °- 2 (3r + 2s a - 2)] , (43) 

where r, s^, s a and £ CT are given by (40). Thus, we 
can calculate the dark energy equation-of-state param- 
eter w\{t) as: 



w A (t) 



PA{t) 
PA(t) 



-16? +8 



k 2 (^ 2s *+i>o 2 i 2s 



(44) 

Finally, expressing (44) in terms of the redshift z we ob- 
tain: 



wa(z) = -- 16^ + 



+8 ■ 



« 2 d0oe 24 ^' n(1+z) -C^o^ 24C " ln(1+2 



K 2 [U<l>y 4U ln(1 + Z) + £ CT Cr ( 2 e- 24 «- ln(l+^)] - 1 



>(45) 



Relation (45) provides w\(z) in terms of the coupling 
parameters £ a and the amplitudes O , <r - Note that 
(45) corresponds to the quintom scenario, and thus ex- 
pressions (40) are embedded in it. Therefore, one cannot 
simply set some parameters to zero in order to obtain the 
simple canonical or simple phantom cases, but he has to 
solve the problem from the beginning with only one field, 
that is the procedure we followed in the previous subsec- 
tions. In the next section we analyze the cosmological 
implications of the quintom model. 



III. COSMOLOGICAL IMPLICATIONS 

In the previous subsections wc have obtained the 
equation-of-state parameter of dark energy w\(z), in 
terms of the coupling parameters £ a and the ampli- 
tudes <j>o, Co- In the present section we investigate the 
cosmological implications for each case. 
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A. Canonical field 

In the case of a simple canonical field, non-minimally 
coupled to gravity, w\(z) is given by relation (16). In 
fig. 1 we depict w\(z) for four different values of the cou- 
pling £0 and for three different values of the combination 
k 2 (/)q. Note that the physical requirement of an expand- 
ing universe, results to an upper limit for £0, namely 
£0 < 1/6, as can be seen in the first relation (12). In 
addition, we mention that in general £0 could be also 
negative, but since it leads to non-physical behavior of 
wa(z) we neglect this case in this subsection. 



N 



2 

< 1 



-1 

-2 
2 

1 



-1 

-2 



1=1/20 



5=1/9 



- - - *V= 1 

■ K** 5 -0.1 

2 

£=1/8 



K \ =10 



conclude that there are some combinations of £0 and K 2 (j)Q 
that lead to a future divergence of wa- Thus, the non- 
minimally coupled canonical field model of holographic 
dark energy predicts a cosmological singularity at a fu- 
ture time tcs, for combinations of £0 and k 2 ^ 2 , that sat- 
isfy: 



(47) 



and since £0 < 1/6, the ^-divergence realization condi- 
tion reads simply: 



>>k 2 4>1 > 1. 



(48) 



Fortunately, this condition leads to a negative effective 
Newton's constant in (4), and thus the corresponding 
parameter combinations must be excluded, leaving the 
model free of a future WA-divergence. In any case, we 
have to mention that in the model at hand the w\- 
divergence at tcs is not accompanied by a divergence 
in the scale factor, in its time-derivative and in the dark 
energy density and pressure. Thus, technically, it does 
not correspond to the Big Rip of the literature [9, 24], 
but rather to some new singularity family. 

For reasons of completeness we present explicitly the 
behavior of wa(z) for k 2 <Pq -C 1, that is for very small 
current value of the scalar field. Specifically, we find that 
the present value wao is: 



FIG. 1: (Color online) wa(z) vs z in the canonical field case, 
for £0 = 1/20, £0 = 1/9, £0 = 1/8, £0 = 1/7, where in each 
case the combination K 2 (f>Q is taken equal to 10, 1, 0.1 respec- 
tively. The divergence of wa(z) is a direct consequence of the 
singularity of (16), and thus the corresponding combinations 
o/^0 and K 2 cj>Q must be excluded. 



As we observe, the value of uua(z) at z = 0, that is 
its current value wao, decreases as £0 increases, while 
its dependence on k 2 </>§ is non-monotonic. However, in 
this simple canonical field case wao is always greater than 
— 1, independently of the values of £0 and n 2 (j> 2 . This was 
expected since this case is well known to be insufficient 
to describe the crossing of the phantom divide wa = — 1 
from above [7]. 

Secondly, we can see that for not so small £0, and for 
k 2 4>q of the order of 1, we obtain a divergence of w\(z). 
This behavior is a clear prediction of relation (16), since 
it possesses a singularity at: 



z s = -l+ (£0K 2 0q) 



(46) 



Therefore, the combinations of £0 and k <f>o that satisfy 
this transcendental equation giving a positive z s , must 
be excluded. However, focusing on the future instead 
of the past, this behavior of Wa has a very important 
cosmological implication. Using directly the form (15), 
which allows us to investigate the future evolution, we 



U 'AOL 2 0=«1 



16^0 > -1, 



(49) 



with the last inequality arising from the upper bound of 
£0 < 1/6. 

Finally, we mention that the model at hand should re- 
ceive additional constraints through the observations of 
the time variation of gravitational constant [27]. In par- 
ticular, differentiating (4) with respect to t and setting 
to = 1 for the present time, we acquire: 



G 
G 



(1-6&)(1-&K 2 $)' 



(50) 



where we have also used (11) and (12). This combination 
must be less than 4% [27]. 



B. Phantom field 

In the case of a phantom field, non-minimally coupled 
to gravity, wa{z) is given by relation (30). In fig. 2 we 
depict wa{z) for four different values of the coupling £ CT 
and for three different values of the combination k 2 <t 2 . 
Note that in this case the physical requirement of an 
expanding universe, results to a lower limit for namely 
— 1/6 < ^cr, as it is implied by the first relation (26). 

As we can see, the value of wao is now a non-monotonic 
function of £ CT and k 2 ctq. Furthermore, we observe that 
for some particular combinations of £ CT and k 2 ct 2 , as a 
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N 



^=1/20 


— k V =1 ° 
— «V=i 




K z c 2 =0.1 
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I 








— k V =1 ° 
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I 





-- K ! O ! =1 




K 2 O 2 =0.1 




^ =1/4 



FIG. 2: (Color online) w A (z) 
for & = 1/20, ^ = 1/8, £ CT 



« z in the 'phantom field case, 
1/6, £<t = 1/4, where in each 



case the combination k o~ is taken equal to 10, 1, 0.1 respec- 
tively. The divergence o/i»a(z) is a direct consequence of the 
singularity of (30), and thus the corresponding combinations 
o/fo- and ft 2 o"o must be excluded. 



consequence of the singularity of (30), there is a diver- 
gence of w\(z) at: 



k 2 ctq, which is not what is expected for a phantom field. 
This behavior is a clear result of the non-minimal cou- 
pling in the holographic dark energy framework. How- 
ever, contrary to the canonical field case where nega- 
tive values of the coupling lead to non-physical behavior 
(ic ao > 1), in this phantom field scenario such a choice 
leads to interesting cosmological implications. In fig. 3 
we depict w\(z) for four different parameter choices with 
negative values of As we observe, negative values 



N 



< 



1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 



■^=-1/6 + , K 2 C7 2 =0.1 

•^=-1/6 + , k 2 ct 2 =0.5 



^=-1/7, k % =0.1 
^=-1/7, k 2 cj 2 =0.5 



1 



(51) 



Thus, the combinations of £ CT and k 2 (Jq that satisfy this 
transcendental equation giving a positive z s , must be ex- 
cluded. Note that in the case of negative condition 
(51) cannot be satisfied and thus this solution sub-class 
is free of a past divergence. 

Similarly to the previous subsection, using the form 
(29), which allows us to investigate the future evolution, 
we conclude that there are some combinations of £ a and 
k 2 ctq that lead to a future divergence of w\. Therefore, 
the non-minimally coupled phantom field model of holo- 
graphic dark energy predicts a a cosmological singularity 
at test for combinations of £ CT and k 2 (Jq that satisfy: 



tcs = (6tK 2 0o) 



> 1. 



(52) 



Note that in the case of negative this condition cannot 
be satisfied and thus this solution sub-class is free of a 
future cosmological singularity. For a positive £ CT the wa- 
divcrgence realization condition is simply 



£ ct k 2 ct 2 < 1, 



(53) 



which does not bring any positivity problems in the ef- 
fective Newton's constant in (18). Thus, the correspond- 
ing parameter combinations cannot be excluded, and the 
model, for a positive £ CT , clearly predicts a future cosmo- 
logical singularity. 

In the case at hand we can see that w\o is always 
greater than —1, independently of the values of <;„ and 



FIG. 3: (Color online) wa(z) vs z in the phantom field case, 
for four different parameter choices with negative values of 
£ CT . The index + in — 1/6 + marks the limit £ CT — > —1/6 from 
above. 



of the coupling produce cosmological behaviors with de- 
creasing w\(z) and w\q very close to —1. Furthermore, 
for £ CT < — 1/7 we obtain a w\q inside the observational 
limits [1-4], although we cannot acquire a clear phantom 
divide crossing. It is known that under specific poten- 
tial choices, a non-minimally coupled phantom scenario 
can achieve the — 1-crossing [28]. It seems that the holo- 
graphic dark energy framework does not allow for such a 
behavior. 

Additionally, taking the limit k 2 <Tq <C 1 we find that: 



WAQ 



«1 



16£ CT > -1, 



(54) 



with the last inequality arising from the lower bound of 
-1/6 

Finally, the present scenario should also receive addi- 
tional constraints through the observations of the time 
variation of gravitational constant [27]. In particular, we 
acquire: 



G 
G 



(1 + 6£,)(1-£^V 2 ) 



(55) 



where we have also used (25) and (26), and thus this 
combination must be less than 4% [27]. 
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C. Quintom model 

In the case of the combined quintom model, that is 
when both the canonical and phantom fields are consid- 
ered to be non-minimally coupled to gravity simultane- 
ously, w\(z) is given by relation (45). In fig. 4 we depict 
wa(z) for four different values of the coupling £^ and for 
three different combinations k 2 (/)q and K 2 a 2 . Note that in 
this case the physical requirement of an expanding uni- 
verse, results to an upper limit for namely £^ < 1/6, 
as it is implied by the first relation (41). The value of 



N 



1- — 



5=1/10 



-k ♦„ =10, 
■ k'4, 2 =1,k(j Z =10 
kV=0.5, tV=».1 



l=M9 



— *x 


= 10, KO ! =1 


---- «\ 






=0.5, k 2 <j„ 2 =0.1 











K ! +„ S =10, K S o/=1 




lc\ ! =1,K ! <s/=10 




k\ 2 =0.5, K Z n Z =0.1 





2 4 



-1 







5=1/6 + 




«\ 2 =10, K Z 0„ 2 =1 




k\ ! =1,k?o/=10 




k\ 2 =0.5, k 2 o 2 =0.1 





2 ^ 4 



FIG. 4: (Color online) wa(z) vs z in the combined quintom 
scenario, for ^ = 1/20, £4, = 1/9, = 1/8, ^ = 1/7, where 
in each case the combinations K, 2 (f>o and k 2 o~q are shown in the 
insets. The divergence of wa(z) is a direct consequence of the 
singularity of (45), and thus the corresponding combinations 
°f £<t>> K^tfro and k 2 o-q must be excluded. 



wao is a monotonic function of As in the previous 
cases, for some particular combinations of £<^, k 2 4>q and 
k 2 (7q, as a consequence of (45), there is a singularity of 
wa(z) at a specific z s . The form of the denominator of 
(45) does not allow for an explicit expression of z s , but 
numerical investigation provides the specific excluded pa- 
rameter values. 

Similarly to the previous subsections, there are some 
combinations of k 2 (/)q and k 2 o"q that lead to a fu- 
ture divergence of wa- Thus, the non-minimally coupled 
quintom model of holographic dark energy predicts a fu- 
ture cosmological singularity for parameter combinations 
that make (44) diverge for tcs > 1, that is in the future. 
We mention that the transcendental form of the denom- 
inator forbids the extraction of an explicit relation for 
tesi but the corresponding values can be provided by 
simple numerical calculations. 

As we observe in fig. 4, wao is greater than — 1, inde- 
pendently of the values of k 2 0q and k 2 (Tq. However, 
for a class of parameter combinations we obtain cosmo- 
logical evolutions in agreement with observations. In fig. 



5 we depict wa{z) for four such combinations of f^, K 2 (j)Q 
and k 2/ Jq. As we can see, we can obtain a decreasing 
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FIG. 5: (Color online) wa(z) vs z in the quintom case, for 
four combinations o/^, K, 2 (f> 2 and k 2 oq shown in the inset. 
The index — in l/6~ marks the limit — ► 1/6 from below. 



form of wa with its current values inside the observa- 
tional limits [1-4]. It is interesting that we cannot ac- 
quire a clear — 1-crossing, which was the basic motive 
of the construction of quintom scenario [10]. It seems 
that the holographic dark energy framework refutes such 
an eventuality. Furthermore, we mention that in the case 
where k 2 0q = k2 °o the effects of the canonical and phan- 
tom fields cancel each other, as expected by relation (45), 
and the dark energy in the model at hand behaves like 
a cosmological constant (dotted curve of fig. 5). Lastly, 
numerical investigations show that the parameter sub- 
space that leads to wao ~ — lj cannot lead to a future 
cosmological singularity, which is also an advantage of 
the model. 

2 ^2^2 



Taking the limit k <PqOq <C 1 we find that: 



^AOU2 2 O .2 <1 



3 - 16^ > 



(56) 



with the last inequality arising from the upper bound of 
£0 < 1/6. Finally, we close this subsection with the ex- 
ternal constraints to the model by the observations of the 
time variation of gravitational constant [27]. In particu- 
lar, we acquire: 



2^2 



(57) " 

where we have also used (40) and (41), and thus this 
combination must be less than 4% [27]. 
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IV. CONCLUSIONS 

In this work we construct various field models of dark 
energy, such as simple canonical and phantom fields, and 
their simultaneous consideration into a combined model 
called quintom. All fields are non-minimally coupled to 
gravity through extra terms in the action, and the inves- 
tigation has been performed in the framework of holo- 
graphic dark energy. In each case we extract wa(z), that 
is the dark energy equation-of-state parameter, as a func- 
tion of the redshift and using as parameters the couplings 
and the amplitudes of the fields, and we analyze it in or- 
der to obtain its cosmological implications. In particular 
we examine the present value wao, the crossing through 
the phantom divide —1, and we extract the conditions 
for a future cosmological singularity. 

For the simple canonical field we find that w\o cannot 
be less than —1, thus this model cannot describe the tran- 
sition through the phantom divide. In addition, we give 
the parameter subspace that has to be excluded since it 
leads to a singular behavior in the past. Furthermore, we 
find that for a specific parameter subspace the universe 
will result in a future cosmological singularity, and we 
extract a specific relation for the time that it is going to 
be realized. Fortunately, the physical requirement for a 
positive effective Newton's constant makes the model free 
of such a singularity. Finally, we give a constraint for the 
model parameters in order for the time variation of the 
gravitational constant to be consistent with observations. 

For the simple phantom field we provide the parame- 
ter subspace that has to be excluded in order to acquire 
a regular evolution in the past. We extract the condi- 



tions and the time of a future cosmological singularity. 
For the case of negative couplings we find a decreasing 
w\ with a current value inside the observational limits, 
in agreement with cosmological observations. The fact 
that i»a lies above the phantom divide is a clear result of 
the non-minimal coupling in the holographic dark energy 
framework. Furthermore, in the single phantom field 
case, the future cosmological singularity cannot be ex- 
cluded. Lastly, we present the constraints to the model 
by the time variation of the gravitational constant. 

For the quintom model, that is the combined case 
of both canonical and phantom fields, we give the 
conditions for physical evolutions, that is without 
divergencies in the past, and we provide the require- 
ments for a future cosmological singularity. We find 
that a clear crossing of the phantom divide cannot be 
obtained, in contrast to what is expected for a quintom 
scenario. It seems that the holographic dark energy 
framework refutes such a behavior. However, we do 
obtain a decreasing w\ with a current value inside the 
observational limits. In addition, these solutions do 
not possess a future cosmological singularity and these 
features make them a good candidate for the description 
of dark energy. Finally we provide the parameter 
constraints in order for the model to be consistent with 
the observed time variation of the gravitational constant. 
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